The box diagram in Yukawa theory 
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We present a light-front calculation of the box diagram in Yukawa theory. The covariant box 
diagram is finite for the case of spin- 1/2 constituents exchanging spin-0 particles. In light-front 
dynamics, however, individual time-ordered diagrams are divergent. We analyze the corresponding 
light-front singularities and show the equivalence between the light-front and covariant results by 
taming the singularities. 



O ' I. INTRODUCTION 

o; 

One of the important theoretical tools to study relativistic bound-state problems is the light-front quantization 
£f) • method which provides one of the three forms [ij of Hamiltonian dynamics and a Fock-state representation at equal 
light-front time r = t + z/c |2|. Distinguished features of light-front dynamics (LFD) may be summarized in two 
fundamental areas of physics, namely the vacuum and symmetry in Minkowski space. Except for zero-modes, the 
vacuum of LFD is much simpler than the vacuum of instant form dynamics (IFD) based on the ordinary equal t 
quantization. The rational energy-momentum dispersion relation in LFD may be the key to understand this highly 
non-trivial property of the vacuum. The number of kinematic generators for the Poincare symmetry is also maximal in 
LFD and the conversion of dynamical generators from the boost to the transverse rotation occurs upon the replacement 
of IFD by LFD . The compact nature of the rotation may be regarded as an advantageous property of the LFD in 
\ handling complicated dynamical operators Q. 

The momentum-space bound-state solutions to a system of relativistic equations tp(x, k±, A) are functions of the 
light-front variables Xi = (fc 4 ° + kf)/(p° +p z ), fe_u, and the particle helicities Aj. They are suitable for the calculation 
of physical observables such as hadron form factors and structure functions. The first step in solving the full set 
Qh. of coupled Fock-state equations on the light-front (LF) is to find a simple analytically tractable equation for the 
valence, lowest-particle-number sector, and to develop a systematic scheme of obtaining the contributions from the 
higher-particle-number sectors with a desired accuracy. A frequently used technique to achieve at least the first step 
is the projection of the manifestly covariant Bethe-Salpeter equation on the hypersurface of LFD jjjj . For instance, the 
the LF ladder approximation can be obtained by projecting the ladder approximation of the Bethc-Salpctcr equation. 
The ladder approximation to the LF two-body bound-state equation involves only up to three-body Fock-states and 
the higher-particle-number sectors beyond the three-body Fock-states cannot be generated in this approximation. 
However, they can be generated if the projection procedure is applied after the covariant Bethe-Salpeter equation is 
iterated once. As an example, the four-body Fock-state of the stretched box diagram in LFD can be generated by 
projecting the covariant box diagram after iterating once the ladder kernel of the Bethe-Salpeter equation. In this 
sense, the procedures of iteration and projection do not commute. 

A three-dimensional reduction of the two-particle Bethe-Salpeter equation has been proposed and the reduction 
of the two-fermion Bethe-Salpeter equation in the framework of LFD has been studied for the 3+1 dimensional Yukawa 
model |7j . We study here a generalization of this model for which the interaction Lagrangian density is given by 

L = + .92*2*202, +ffl2(*l*2CT + h.C.), (1) 

where the fermions correspond to the fields *i,2 with rest masses M, m and the exchanged bosons to the fields 0i, 
02, and a with masses fix, /i2, and respectively. A physical system where this interaction would be applicable is 
the coupled NN - N*N* system. 

The LF treatment yields three-dimensional quantities for the transition matrix and the bound-state wave function. 
Since the kernel of the LF Tamm-Dancoff reduced Bethe-Salpeter equation for the vertex function for the one- 
boson- exchange interaction in ladder approximation has a divergence problem, the introduction of a counterterm to 
renormalize the integral equation was proposed In Ref. 0, the authors were concerned with the origin of the 
perturbative counterterms of the LF ladder Bethe-Salpeter equation for the Yukawa model and showed that the kernel 
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of the auxiliary integral equation, expanded up to the fourth power of the coupling constant, g , naturally yielded 
the box counterterm [t| and a well defined finite part. As the authors explained, this is because the perturbative 
expansion of the LF scattering amplitude in powers of the coupling constant, obtained from the LF T-matrix equation 
with the kernel calculated up to the same order, necessarily reproduces the perturbative covariant ladder scattering 
amplitude at that order in g. An important necessary inclusion, in addition to the instantaneous contributions, is the 
higher Fock components by going beyond the Tamm-Dancoff approximation. For example, in Ref. the divergence 
could not be removed if the four-body Fock-component represented by the simultaneous propagation of two as and 
two fermions occuring between the creation and annihilation of the bosons was not included. Even higher Fock-state 
contributions such as the simultaneous propagation of three as and two fermions should be necessary if terms of order 
g 6 are included in the kernel. Such a requirement of including higher Fock-state contributions to remove the relevant 
divergence in LFD would persist if higher and higher Fock sectors would be included. Thus, it is not clear how the 
reduction program would work in practice even if the expansion in powers of the coupling constant allows a definit 
number of boson exchanges. 

Moreover, LF perturbation theory shows singularities which have nothing to do with the truncation of Fock space. 
Many pitfalls and treacherous points exist in LFD, which is full of surprises. In particular, we rece ntly discussed 
the arc contribution and the point singularity in the contour integration of the LF energy variable llOl . Also, we 
showed an explicit example of the anomaly associated with the quantum field theoretic infinities that make 
the prediction of physical quantities different between the manifestly covariant approach and LFD even after the 
amplitude is renormalized, unless the anomaly-free condition is imposed. Our vector anomaly analysis 11] provided 
a bottom-up fitness test of the Standard Model and a model-independent proof of the zero-mode contribution even in 
the good (or plus) current matrix element (helicity zero-to-zero amplitude). Thus, it is significant to further analyze 
the LFD singularities. 

In this work, we use the generalized Yukawa model given by Eq. Q as a testbed for different ways to remove 
LF singularities and discuss another type of singularity which is significant to the bound state problem in LFD. In 
LFD one would like to use the LF time-ordered one-boson-exchange amplitudes as driving terms in the bound-state 
equation, similarly to the covariant one-boson-exchange amplitudes which play the same role in the ladder Bethe- 
Salpeter equation. Here a difficulty arises that can be simply formulated in perturbation theory: The manifestly 
covariant box diagram in the Yukawa model is finite, whereas the corresponding LF diagrams are divergent. In this 
type of diagrams, the residue calculus is correct and no arc contributions occur in the LF energy contour integration. 
Moreover, the LF divergences cancel [l2| since this type of singularities corresponds to finite integrals over the LF 
energy variable. However, in order to learn how to exploit the property of cancelling divergences, we analyzed in 
detail the box diagram. We present both the manifestly covariant calculation and the LF calculation to verify the 
equivalence between the two. Using the integrals over the Feynman parameters in the manifestly covariant calculation, 
we expand the on-shcll scattering amplitude in terms of form factors that are functions of the Mandelstam variables. 

In order to be able to calculate the LF amplitudes, we need to introduce regularization. Two methods were used, 
namely Pauli-Villars (PV) regularization where one PV boson was introduced, and dimensional regularization in the 
transverse- momentum integrals (DR2). We compare these two methods, because we have seen before 1 1 1 1 that in 
some cases the finite parts of the regulated amplitudes may differ. Our LF calculation reveals the zero-mode in the 
stretched box depending on the kinematics. We discuss a few different prescriptions to compute the form factors. We 
have found that they are prescription independent. 

In Section II, we briefly review the scalar box diagram and introduce the variables that are used in the rest of the 
paper. In Section III, our manifestly covariant calculation is presented for the scattering of two spin-1/2 particles. 
The relation of the matrix elements in spin space to the invariant form factors is given here too. We show how, for 
a given kinematics, the form factors can be extracted from the matrix elements. In Section IV, our LF calculation 
is presented and the zero-mode contribution from the stretched box is also discussed. In Section V, we sketch the 
two regularization methods we use in the LF case. Section VI contains a discussion of the numerical results and our 
Conclusions follow in Section VII. In the Appendices A and B, details of the helicity spinors and the integrals used 
in DR2 are given, respectively. 



II. ANALYSIS OF THE SCALAR DIAGRAM 

In Refs. ^3 an d 01 > a general framework is given for the calculation of the scalar box diagram. Although these 
authors were chiefly interested in the analytic properties of amplitudes in strong-interaction theory, we may use their 
methods to determine the values of the Mandelstam variables for which the box diagram is non-singular. 

The scalar box diagram is defined by the integral 

f 4 , f , S(l — ot\ ■ ■ ■ — a.4) . s 

1= a k dai . . . don- — ^ ^- ■ r - 2 jTfj. (2) 

J J [ai(q(-m() + h a 4 (g| - m|)] 4 




p 2 ' p,+ p 2 - *■ p 2 X 

FIG. 1: The covariant box diagram with momenta defined. 



Of course, Feynman parametrization was used here. The momenta qi on the internal lines are (see Fig. ^) 

qi = k-pi, q 2 = k, q 3 = k - p[, q 4 = pi + p 2 - k. (3) 

The corresponding masses are 



mi = TO3 = /i, TO-2 = TO4 = m. 

The function in the denominator can be written as k' 2 + M 2 ov if the shift 

k — > k = k — [a\ + oti)pi — a^p2 



(4) 



(5) 



is performed. 

Next, we make the restriction that the external momenta are on shell for external particles with mass M. Owing to 
four-momentum conservation we can express the momentum dependence in terms of the three Mandelstam variables 
s, t, and u, defined as usual: 



s = (P1+P2) 2 = (P'i+P'2) 2 , 
t = (ti-Pi) 2 = (P' 2 -P2) 2 , 

U = (Pi -P2) 2 = (P[ -P2?- 



(6) 



Upon substitution of the on-shell relations for the external momenta and the Mandelstam variables s and t we find 



M c 2 ov = {a>2 + a 4 )m 2 + (a 1 +a 3 )fi 2 

— {ct\ + Qf3)(a2 + a.i)M 2 — a.2aiS — a\a 3 t. 



(7) 



We consider in this work a situation that is similar to the bound-state case in the sense that the corresponding 
scattering amplitude is real, namely we consider masses m and M, related to the fermions in the loop and the 
external fermions respectively, such that M < to, which gives us a window AM 2 < s < 4to 2 where the contribution 
from the ^ ^-intermediate state to the scattering amplitude is real. While working in this window, we can avoid the 
complications caused by the unitarity cuts, when we calculate the diagram for the process ^i^i — + ^2^2 —> ^i^i. 
In that case the analysis of Refs. |13L FLU ] shows that the minimum value of the denominator function defined above 
is greater than zero in the domain of integration. 

III. COVARIANT CALCULATION OF THE YUKAWA BOX 

In this section we discuss the general formalism for the calculation of the box diagram. First, we sketch the 
calculation of the amplitude in spin space and next we connect these matrix elements to invariant form factors. 



A. Amplitude 

We consider the scattering of two spin-1/2 paricles with mass m, that exchange scalar particles of mass /i. The 
amplitude can be written as 



Tfi = u(p x , s[)u(p' 2 , s' 2 )Mu(p 1 ,s 1 )u(p 2 , s 2 ), 



(8) 



4 



where M. is a matrix in spin space and depends on the invariants that can be built from the momenta px, p 2 , p'i, and 

The spin matrix Ai, defined in Eq. (JHJ, can be obtained using the standard Feynman rules. Then one finds 

f d 4 fc (ft + m) (g> {fa + p* 2 - Jf + m) 
J (2tt) 4 (k 2 - m 2 + it)((k - Pl ) 2 - ^ 2 + it) 
1 



({pi +P2- k) 2 - m 2 + it)((k - p[) - [i 2 + it) ' 

where the notation ® is introduced to distinguish the internal lines connecting pi and p[ for i — 1 and 2. This 
expression can be rewritten using Feynman parameters as 

M^j TAai - Aai J ^ + m ^'*^-* + m \ (10) 

where 

D = a x ((k - Pl ) 2 - Ai 2 )) + a 2 {k 2 - to 2 ) 

+a 3 ((k - pi) 2 - /i 2 ) + 04((pi + p 2 - k f - m 2 ) + it. (11) 

The domain of integration T is the interior of a regular tetrahedron given by T — > 0, a\ + a 2 + a 3 + a 4 = 1}. 

After a Wick rotation it is clear that the amplitude is given by a convergent integral. Consequently, we are allowed 
to shift the integration variable. The shift is of course the same as in the scalar case, Eq. (0. After the shift, which 
produces a denominator that is an even function of the integration variable fc', we may invoke symmetry to prove 
that the terms in the numerator proportional to k' vanish upon integration and that the term k'^k' v is equivalent to 
a term g^ v fc' 2 /4. 

The numerator is changed to (we drop the terms linear in k') 

-7(1) • k! 7(2) • k' + {7(1) • [(1 - as - a,i)pi + a 2 p 2 + a 3 p[] + m} {7(2) • [(a 3 + a 4 )pi + (1 - a 2 )p2 - a>3p[] + m} . 

The gamma-matrices are associated with two different particles: 7(1) is associated with the internal line connecting 
Pi and p\. 

We need the following integrals 

D = 6 / da 
Jt 

D a — 6 da 



T 



D aiaj = 6 I da 

T 



T 



d 4 k 


1 


(2tt) 4 (k 2 


- M 2 ) 4 ' 


d 4 k 


on 


(2tt) 4 (k 2 


— M 2 ) 4 ' 


d 4 k 




(2tt) 4 (k 2 


— M 2 ) 4 ' 


d 4 k 


k 2 


(2tt) 4 (k 2 


— M 2 ) 4 



D 2 = 6 da -—j—, -= -. . (13) 



Upon performing the integral over the momentum and using a Wick rotation, we find 

B n , = \ 9 / da- ' 



(An) 2 J T M 4 



COY 



D, y . n — -r. — err / da- 1 



{Ait) 2 J t M 4 



COV 



^wrL^ik- (14) 

Owing to the symmetries of the diagram the mass function is symmetric under the transpositions a\ «-> a 3 and 
a 2 <-> a.4. Using this symmetry of the denominator, we find the identities 



D 


= 2{D ai 


+ D a2 ) 7 


D ai 


— Da 3 j 




D a{ 


= D <4> 


D al =D al , 









(15) 
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Using the fact that the alpha's add up to 1, we can derive two more relations, namely 

-D a 2 — F) a2 2_D aiQ2 F) CX2CXA . (1^0 
Using the symmetries just discussed, we find that only D ai , D a2 , D aia2 , 73 QlQ3 , 7J> Q2Ct4 , and Di are independent. 

B. Spin Structure 

We now gather the different pieces of the amplitude. We write 

M = -7(1) ■ 7(2) \D 2 + m 7 (l) • D{1) + m 7 (2) • D(2) + j(l)^(2) v D(12)^ + m 2 D , (17) 
where we used the vectors 

D(l) = | [PD + qD ai - q'D ai ] , 
D{2) = | [PD Q - qD ai + q'D ai ] , 

= PD -D(1), (18) 

and the tensor 

D(12)»» = ~P»P»[D ai +2D aia2 +2D a2ai } 

+ \ W + tptf] [~D ai + 2D aia2 + D aia3 ] 

-l P ^D ai - ~P^D ai - ~q^q'^D aia3 . (19) 

Here we have used the vectors P = p\ + p%, q = p± — P2 and q' = p[ — p' 2 to get symmetric expressions. We use the 
notation p^q"} = v ^ q v + q^ v v an d p&q u ] = p^q" - q^p v . 

In view of the symmetries of the spin matrix elements of the identity and the 7-matrices (see Appendix^, we find 
that the matrix T has the following structure 





Tia 


T 13 


T u 


T 21 


T~22 




T~24 


~ 2 24 




2 22 




2 14 


~ 2 13 


~ 2 12 





2 I ST* ry* rf*. /t-* I ■ J 



Here we use the following numbering of the two-fermion spin states 

|1) = ITT),|2) = |U) ) |3) = UT),|4) = |U). (21) 

Clearly, there are eight independent complex matrix elements which correspond to sixteen independent real numbers. 

Upon taking matrix elements between spinors and using the Dirac equation to simplify some matrix elements, we 
find the following structure 

T = 1 F 1 + 2 F 2 + 3 F 3 + 4 F 4 , 

1 = 7(1)® 7(2), Fi = 2mM {2D ai + D a2 ) + M 2 (2D ai + D a2Ct4 ) + 2m 2 (D ai + D a2 ) , 

2 = [7(1)8 (pi-r(2)) + (p2-r(l))® 1(2)], F a =mD aa +M(2D ataa +D aacu ), 

3 = (p 2 -r(l))®(p 1 -r(2)), F a = D a2Cti , 

4 = T(l)-r(2), F 4 = -73 2 /4. (22) 

Apparently there are four independent form factors. The reason for this small number is that there is a high degree of 
symmetry in the mass function M 2 ov . If this symmetry would be broken by e.g. differences in masses of the particles, 
more form factors would occur. Remarkably, the integral D aia3 does not occur. 

Using this representation, one can easily derive linear relations between the spin matrix elements if the matrix 
elements of the spin operators O 1 . . . O 4 are known. 
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C. Extraction of form factors 



In a manifestly covariant calculation, one does not need the spin-matrix elements Tfi to extract the form factors 
as Eq. IL'l't shows. However, in a LF calculation, which breaks manifest covariance, there does not exist a relation 
like Eq. Q22)). Thus, one must try to solve the linear relation between the Tfi and the Fj in a specific kinematics. We 
work in the center of mass system (CMS). Of course, other kinematics will also enable us to extract the form factors, 
but the CMS, which is detailed below, seems to be the simplest one. 

For any kinematics, we may write 

■J*. »:Mv ( 23 ) 

3 

The extraction procedure consists in choosing a set of spin labels {/z,«z}, I — 1, ... ,4, such that the square matrix 
J f L 4 , j = 1, . . .4, I — 1, . . . , 4, is nonsingular. An obvious choice is fixing the row label / to 1 or 2 and let the column 
index i run from 1 to 4. This choice was made in this work and we have found that it gives unambiguous results as 
long as the sine of the scattering angle 9 and the three-momentum p of the scattered particles do not vanish. 

As an illustration of the problems one may encounter otherwise, we show the structure of the amplitude matrix in 
forward [6 = 0) and backward (9 = it) kinematics. In forward kinematics, the amplitude T reduces to a diagonal 
form 



T n 














Tn 














T n 

















T(9 = 0)=\ Q » Tn Q . (24) 
V Tn J 

Clearly, there is only one independent amplitude in forward kinematics and one obviously cannot extract four form 
factors at 9 = 0. 

In backward kinematics, the amplitude T is a slightly more complicated and has the following structure 



Tn 














Tn 


T23 








T23 


Tn 














Tn 



T V = «)=\ oZTn ■ ^ 
V Tn J 

Now two independent amplitudes occur, which also does not allow to extract the four form factors. 

In the LF case, we may use the same matrices O 3 ^ to connect the LF amplitude corresponding to diagram d 

to an LF 'form factor' Fj . We use here the terminology of form factors, although the quantities Ff are not invariant 
objects. We shall, however, show that upon regularization of the LF amplitudes, the corresponding LF form factors 
add up to the invariant ones, viz. 

Y. F t = F r ( 26 ) 



D. Kinematics 



In the box diagram that depends on three independent momenta, a simplification cannot be achieved by choosing 
e.g. the Breit frame, which appeared to be so helpful in, for instance, the triangle diagram. Therefore, we choose just 
the CMS. A slight simplification can be achieved if one limits the incoming and outgoing momenta to the xz-plane. 
Here we follow Erkelenz 18]. In particular, we choose 

Pi = pe z , p{ = p(sm9 cos <fi, sin 9 sin0, cos 9), p 2 = —Pi, p% = —p{- (27) 

In order to construct the correct helicity spinors in initial and final states, it is important to define the polar and 
azimuthal angles correctly. They are 

Pi : 0i = 0, <f>x = 0; p 2 : 9 2 = vr, 2 = tt, p\ : 9[ = 9, </)[ = 0; p 2 : 6' 2 = tt - 8, 4>' 2 = w. (28) 

The corresponding spinors can be found in Appendix ^ 

Here, we want to point out that the zero mode we mentioned in the Introduction is found in this kinematics when 
p^ = p'^~ , which occurs for p — or 9 — 0. If we had chosen instead of the ccz-plane the xy-plane as the scattering 
plane, the zero mode would have occurred for any values of the scattering angle and the three-momentum. 
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mi 



(in 



FIG. 2: The splitting of the covariant box into light-front time-ordered diagrams. 



IV. LIGHT-FRONT CALCULATION 

The invariant amplitude Eq. I jlOj l can be rewritten in LF coordinates 

fd 2 fc ± fdk+ f dfc- 1 {ty + m)®{-$ + j> l +j> 2 + m) 

™ J (2tt) 2 J 2ir J 2tt $ (k- - ffr)(fc- - ff a )(fe- - tf 3 )(fc- - ff 4 ) ' 1 j 

where the phase-space factor $ is given by 

$ = Wk + {k + - pf)(k + - p[ + )(k+ - pf - P +). (30) 

The 'Hamiltonians' Hi are defined as 

(k 1 - - pi) 2 + p, 2 - it 



Hi = Pi 
(k- 



2{k+-p+) 



H 



2 



H 3 = p\ 



i — 



2k+ 

(k 1 - -p'^) 2 + p 2 -it 



2(k+-p[ + ) 
(k^ -pj-pjr) 2 

2(k+-p+-p+) 



Ha = Pi +P 2 H ^77Z + — ~R • ( 31 ) 



We choose p± > p[ + for definitness. In the opposite case, some details of the calculations would change, but the 
general picture, and in particular our conclusions, would remain the same. 

In order to get the light-front time-ordered (LFTO) diagrams, we integrate over k~ first. The positions of the poles 
define five regions in k + ', each of them with a different number of poles in the upper and lower plane. They are 

1. k+ < : Imifi > 0, lmH 2 > 0, lmH 3 > 0, Imif 4 > 

2. < k+ < p[ + : lmH 1 > 0, Imif 2 < 0, ImiJ 3 > 0, lmH 4 > 

3. p[ + <k+ <p+ : ImHi > 0, lmH 2 < 0, Imi? 3 < 0, Imif 4 > 

4. p\ < k+ < p+ + p+ : JmHx < 0, ImiJ 2 < 0, lmH 3 < 0, ImH 4 > 

5. k+ > p+ + p+ : Im-ffi < 0, ImiJ 2 < 0, lmH 3 < 0, lmH 4 < 0. 

In the regions 1 and 5, all poles are in one half of the complex plane so that the integral over k~ vanishes. In regions 
2 and 4, one of the poles is in one half plane, while the other poles are in the other half plane. We close the contour 
in the half plane with one pole and perform the integration. In region 3, there are two poles in either half of the 
complex plane so that for this calculation two residues have to be included. 

Since two of the internal particles are fermions, the instantaneous parts have to be taken into account. We obtain 
the LFTO diagrams of Fig. |3 when we perform the described analysis. Diagrams (a) and (b) correspond to region 2, 
diagrams (c), (d), (e), (f), and (g ) correspond to region 3 and diagrams (h ) an d (i) correspond to region 4. We can 
decrease the number of diagrams we have to use with the blink mechanism |l9j [. which we will do for each region. 
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Region 2 

We can combine the propagating diagram with the instantaneous diagram in region 2 to obtain a diagram with a 
blink, as shown in Fig- El We will call this diagram with the blink diagram 1. 



Sj t- 




FIG. 3: The sum of the two diagrams of region 2 is a diagram with a blink. 



The corresponding amplitude of diagram 1 is 

fd 2 k^ dfc+ -i [# on + to] ® [ 7 „ K + - Kn) + m] 

where the momenta with subscript 'on' are four-momenta with the minus-component calculated from the on-mass-shell 
condition, i.e., k~ n = ((fc- 1 ) 2 + m 2 )/2k + . The energy denominators are 

; _ (k^r+m 2 (p^-k^r+^ 

JJ\ = tl-i — Hi = p x - 



D 2 = H A - Hi = Pl + p 2 



2fc+ 2 (p[ + k+) 

(k^) 2 + rn 2 (p^ +p^ -k^+m 2 
2k+ 2(pf +pt - k+) 



_ (k^) 2 +m 2 (pt-V^+j* 

D 3 =H 3 -Hi= Pl -^^ 2(pi-k + ) ' (33) 

The phase-space factor is 

$ = I6k+(k+ - p+){k+ - p' 1 + ){k+ - p+ - p+). (34) 



Region 3 

In this region two poles contribute, which are not directly coupled to diagrams, see Ref. 0| f° r a discussion. Each 
of the diagrams is a linear combination of the corresponding residues. The diagrams contributing in this region are 
(c), (d), (e), (f), and (g) in Fig. |2| The open diamond, diagram (c), does not contain any blinks. We will call this 




FIG. 4: The sum of the instantaneous diagrams of region 3 and the stretched box is the stretched box with two blinks. 



diagram 2. Its amplitude is given by 



'd 2 ^ ft dk+ i [$ on + to] g Mpi +P2- ky on + to] 
M2 ' W? i+ ~W ¥ D^Dl • (35) 
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The instantaneous diagrams (e), (f), and (g) can be combined with the stretched box diagram (d), to create a stretched 
box with two blinks (see Fig. 0J|. The corresponding amplitude is 

d 2 £^ ft dk+ i YlM - (pi - + m] ® [ 7l ,(p^ - (fc - K)ob) + "»] 

where 



^ = iwlt« ' (36) 



(PT + - fc^) 2 + ™ 2 ~ Pi' X ) 2 + M 2 



D[ = H 4 -H 2 = p' 2 

D 2 = H 3 - H 2 = p 2 -p 2 7 ^ TZ — — . (37) 



2(p++p+-k+) 2(k+-p'+) 

(fc x - p{^) 2 + /i 2 _ (p^-fc^T^ 2 
2(fe+-p'+) 2(p+-fc+) 



We will call the stretched box with blinks diagram 3. This diagram has a subtlety, i.e., it contributes even when the 
integration region goes to zero; this is called a zero mode. In Sec. llVCJl we will elaborate on this subject. 

Region 4 

In this region, only the pole corresponding to H4 contributes to the amplitude so that we now have diagrams (h) 
and (i) in Fig. El Again we can sum these two diagrams to obtain the diagram with the blink (see Fig. [21 ■ We will 
refer to this diagram as diagram 4. The corresponding amplitude is 

fd 2 ^ M+Pj d fc+ -i [yM+P% - (pi +P2- fcjftj + m] <E> [ lv { Pl + P2 - kf on + m] 

Mi -Jwf] Pt • (38) 

The energy denominator different from those of regions 2 and 3 is 

n « „ rr - (pj+pj-k^+m 2 (k^-pjf+l? 

D 3 = Hi- Hi=p 2 , + — — +— . (39) 

2{Pi+P 2 -k+) 2(k+-p{) 

A. Amplitude 

When we take the matrix elements of the four Mi derived above, between spinors, we end up with the contribution 
of the diagram to the amplitude T . All four diagrams have the same structure so that we can parametrize the 
contribution of each diagram to the total amplitude T in the following way (d labeling the four LF contributions) 

T d = u(p' 1 ,s' 1 )u(p' 2 ,s 2 )M d u{p 1 ,s 1 )u{p 2 ,s 2 ) 

= c(d)^r(i) M <8 r(2)„ + c(d)^r(i) M ® 7(2) + c{d) v iv i{i) ® r(2)„ 

+C(d) u I(l) <g> 1(2). (40) 



FIG. 5: The sum of the instantaneous and propagating diagrams in region 4 is equal to a diagram with a blink. 
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The tensor C(d)p r , the vectors C(d)p 7 and C(d)j r and the constant C{d)u are of course different for each diagram. 
Taking for example diagram 1, we obtain 



C(l 



J (2tt) 2 J 2tt $ Di£> 2 £> 3 
d 2 /^ /• p '^dfc+-i fc£,m 



C(1) " I (2tt) 2 ,/ 2tt $ D °DM 



cm 



IT 



i + 



d 2 fc^ f p 1 dfc+ -i m (K + - fc£ 



(27t) 2 J 2tt $ Di£> 2 £> 3 



The expressions for the other diagrams are similar. 

In principle, we now only have to calculate C(d)p r , C(d)p 7 , C(d)j r and C(d)u for each diagram. Using Eq. (|40|) 
we obtain the corresponding T d and then use Eq. I|23[) to extract the form factors. We end up with the contributions 
of each diagram to the form factors. 

B. Rationalization 

Because the numerator and denominator contain terms like l/(p'^~) 2 , it seems there are singularities for the k + 
integration. But this is not the case and to show this, we rationalize the denominators. Let us, for example, take the 
first diagram and define 



di = 2fc+ (pi T - fc+J D u 
d 2 = 2k+ (p+ + p+ - k+) D a , 

d 3 = 2k+ (p+ - k+) D 3 . (42) 

These d^'s are rationalized and we can proceed to obtain rationalized denominators for each diagram. We now use 
for diagram 1 

I— L_ = M)!_L_ (43) 

$DiD 2 D 3 2 did 2 d 3 V ; 

to rewrite all the integrals in terms of the d^'s. Then, it becomes clear that the /c + -integration is not singular. 

C. Zero modes in the box diagram 

Domain 3, p l+ < k + < p + , corresponds to two LFTO diagrams; the open diamond and the stretched box. If we 
combine the LF-propagating and instantaneous parts to blinks [19|, the longitudinal singularities cancel out. It is easy 
to show that the open diamond is obtained by integrating a function that is finite in domain 3. Thus, its contribution 
vanishes if this domain shrinks to a point. The stretched box, however, contains a zero mode. 

To show this, we look at the integrand for the stretched box which is given by Eqs. (|33() and (|34() . Upon rational- 
ization, we find that the integrand is proportional to Tz given by 



l (p+-k+)(k+- p'+)N 3 
2 d'^dz 



where iV 3 is the numerator for A4 3 . The factors (pf — k + ) and (fc + — py) are just the ones needed to cancel the 
longitudinal singularities in 7V 3 . Indeed we find 

(p+ - k+) [ 7m K - (pi - *)£,) +m] = - fc x ) 2 + M 2 ]7 + (l) + (Pt - k + ) x finite terms, 

_ p '+) [ lv{p ' 2 v -(k- p'^J + m ] = ~[{k x - p'^f + fi 2 h + (2) + (k+ - p'+) x finite terms. (45) 
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Thus, we see that upon taking the limit p± — > p± the numerator reduces to 

= i 7 +(l) 7 +(2)[(pt - +M 2 ][(fe X -P1 X ) 2 +M 2 ]- (46) 

Next we take the same limit in the denominators to find the final result for the zero mode integrand 

^ zm = _ 7 + (l)7 + (2) 1 m 

8k+(p+ +p+ - k+) (k+ -p'+Wt - k^) 2 + M 2 ] + (pt - fc+))P ± -P'^) 2 + M 2 ] ' 

We see that the zero mode occurs in the 7 + (l)7 + (2) tensor element only. Moreover, the logarithmic divergence is 
clearly visible: for lAr 1 ! — > oo, !F zm behaves like 1/k ±2 . Finally, the factors p+ — k + and k + — p'-^ in the denominator are 
just the factors that cause the zero mode to survive. The integral over k + can be done by making the transformation 

k+ =p'+ +x(p+ -p'+). (48) 

Then the zero-mode amplitude is obtained after taking the limit p± — ► p± . It is 

^=-3^^ r dx /w , l_ . (49) 

8p+p+ (2tt)3 J J x {p± - k^f + (1 - x)(k x - p{^) 2 + /i 2 

The divergent integral over the transverse momentum can be regulated using DR2 or PV. 



V. REGULARIZATION METHODS 



Now that we have the expressions for the four diagrams, we can see that there are ultraviolet logarithmic divergences 
in the ^-integration. These divergences have been explored by Van Iersel |2l]]. The four diagrams all have a term 
proportional to (fc -1 ) 4 in the numerator and terms proportional to (fc^) 6 in the denominator, making the integral 
logarithmically divergent. Van Iersel has shown that the divergences cancel when one sums all the diagrams, but the 
question whether the summation of finite parts gives the covariant results remained open. We discuss two possible 
regularization schemes, two-dimensional dimensional regularization (DR2) and Pauli-Villars (PV) regularization. 



A. Dimensional regularization 

In the LF case, we only need to regulate the transverse directions, so that we compute the diagram as an analytic 
function of the dimensionality D of these directions. The final quantity should then have a well-defined limit as 
D — > 2; we call this form of dimensional regularization two-dimensional dimensional regularization (DR2). 

The divergent parts correspond to terms in the numerator proportional to fc 4 . In order to use DR2, we have to 
rewrite the denominator in the form D = ((fc^) 2 + (M ± ) 2 ) 3 . We will obtain this form by introducing Feynman- 
parameters followed by a shift of the integration variable. 

The denominators of all the diagrams have the same structure, D = dic^cfa- We introduce Feynman parameters to 
change the numerator into 

, , , = I da 1 ---a 3 - -3. (50) 

«1«2«3 Jo \QLld\ + Q!2G?2 + Q!3G?3) 

To obtain the desired form, we expand the denominator in terms of k 1 - as follows 

D d = (a d (k^ 2 + bj ■ k 1 - + c}j , (51) 

which can be rewritten in the form 
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where for a diagram labeled d we write 



{Mi 



2a d 

2 



Cd 

ad 



(53) 



We shift the integration from k 1 - to k and also express the numerator in terms of k . Then, using the usual 
symmetry argument, we see that only even powers of k' 1 - in the numerator contribute and terms in the numerator of 

The final expressions after the shift are rather 



the form (jz ,J - ■ a 1 - J (k' ± ■ b ±S j can be replaced by i 



■ Ft 



complex and thus we will not give them here. 

Finally we conclude that we can write C(G?)pp, C(d)£j, C{d) v IT and C(d) II in terms of the following integrals 



r - 



I = 



p 



d 2 k 




1 






) 2 + (M-L) 2 ) 3 


d 2 k 






(2tt; 




) 2 + (M- 1 ) 2 ) 3 


d 2 k 


J_ 


(^) 4 


(2tt; 


2 ((k x 


) 2 + (Af^) 2 ) 3 



(54) 



After the /c^ integration is done, we also have to integrate over k + and the Feynman parameters. 
Since the integral if is divergent, we regulate it. Using DR2 we write the integral as 



,2-D 



f'-s 



d°k A 



{k^f 



(2tt) 



((fc-L) 2 + (M J 



\3 ' 



(55) 



where D = 2 — 2e and a scaling mass /i s is introduced for dimensional reasons. This integral can be done as shown 
in Appendix [BJ Finally we end up with the following expression for this integral 



-2 1 

T 2 — 

73 _ in 



1 

7' 

e 



- + In (4vr M 2 ) 



In M 



(56) 



The C(d)pp, C(d)pj-, C{d) u IT and C(d)u are linear combinations of these integrals. When the - — 7 — | + In (47r^ 2 ) 
terms of the four diagrams are summed, they cancel out, so we will not take them into account. The cancellation is 
proved below. 

The divergent terms are the terms in the numerator which are proportional to (fc'^) 4 . Only the C(d)^r contain 
such terms. We will now show that the divergent parts of these terms cancel. The divergent parts of the four C(d)pp~'s 
for d = 1, . . . , 4 can be written as 



Cd 



1 

- - 7 " 



- + In (4^ 2 



(57) 



with the quantities Cd given by 



c 2 



I 

47T 



p 'l dk+ 







I 

47T 



2tt 

p * dk+ 
,+ ~2tt~ 



da 



A 



da ■ 



1 

1 
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where the quantities ad are defined in Eq. (|51|l . The notation J^da is used for the integration over the Feynman 
parameters, which lie in a triangle in a-space. 

If we substitute the values found for ad, perform the integration over the Feynman parameters and sum the four 
integrals, we get 



i 1 

4^8 



1 

- - 1 
t 

p " dk+ 



In (4tth 2 s 
1 



p '* dk^ 



1 



2tt - 

p i dk+ 



-p'i (pt + pi) pi 



,+ 2n -p'+ (p+ + p+) p+ J p , t 2vr -p'+ (p+ - p'+) p+ 
Pt+vt dk+ 1 



2tt 



-p' 2 + (pi + pi) pi 



(59) 



Performing the /c + -integration, we find that the factor multiplying [l/e — 7 — 3/2 + In (47r/i 2 )] vanishes, so the diver- 
gent terms of the diagrams cancel out. 



B. Pauli-Villars regularization 

In the previous section, we used DR 2 to regulate the divergent integrals of the four LFTO diagrams. This was 
only necessary for the LF calculation, since there are no singularities in the manifestly covariant case. There are of 
course other ways to regulate the divergent terms, for example Pauli-Villars regularization, which is the topic of this 
subsection. 

In the case of Pauli-Villars regularization, we change the meson propagator in the following way 



p 2 — fi 2 + it p 2 — (i 2 + it p 2 — A 2 + it 

which means we are adding an extra particle to the theory. All the integrands are unaffected by this change at low p 
(since A is large), but they are cut off smoothly when p > A. 

We need only one Pauli-Villars particle because we are only regulating the box diagram. If we want to regulate 
more diagrams, we may need more Pauli-Villars particles (see for example Brodsky et al. |2fj|). 



VI. RESULTS 

We evaluated the amplitudes using the two regularization schemes given the in Sect.0 We found that the covariant 
calculation gives the same results as the LF calculation, i.e., they are equivalent. In the calculations, we have used 
the following values for the masses: M = 0.94, m = 1.44 and /i = 0.14 corresponding to the masses of the hadrons 
N, Roper resonance, and 7r. These values satisfy the stability conditions of Sec. [H] 

In all our numerical calculations, we have used Gauss-Legendre quadratures to perform the integrals over the a's 
and k + as well as an adaptation of Gaussian quadrature for the integrals over k . We have checked the convergence 
of all numerical calculations. 



A. Manifestly covariant calculation 

The matrix elements and the form factors depend on two independent quantities, namely the Mandelstam variable 
s (which is directly related to the total energy) and the scattering angle 9. 

First we give the results for both kinematics possible for 7n in 1+1 dimensions. These results are shown in 
Fig. |SJ The lines intersect at the lower threshold s = 3.53. This has to be the case, since at this threshold forward 
scattering cannot be distinguished from backward scattering. For comparison, we give in Fig. [3 the matrix elements 
of T±i in forward and backward scattering for the 3+1-dimensional case. The behaviour is almost the same as the 
1+1-dimensional case, the main difference being the scale. 

Of course, we also have in the 3+1 dimensional case the 723-element in backward scattering, but as its behaviour 
is very similar to the behaviour of Tn, we do not show it. 

In Fig. [SJ we present our results for the s- and ^-dependence of the four form factors. 
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FIG. 6: The s-dependence of T\\ in forward and backward kinematics in 1+1 dimensions. Covariant calculation. 
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FIG. 7: The s-dependence of Hi in forward and backward kinematics in 3+1 dimensions. Covariant calculation. 



B. Light-front calculation using dimensional regularization 



We have calculated the s-dependence of the contributions of the four diagrams to T\\ in 1+1 and 3+1 dimensions 
both in forward and in backward kinematics. The results are given in Figs. - El Note that in these cases the 
contributions to the form factors cannot be calculated, as the matrix Oj i of Sec. IIII CI is singular. 

In Figs. I§1 and 1111 the zero mode is clearly visible. In 1+1 as well as in 3+1-dimensions, diagram 1 dominates 
at high s in forward kinematics. In backward kinematics, it is diagram 2 which contributes the most. We also see 
differences between 1+1 and 3+1 dimensions. The dominant diagrams are the same, but otherwise the behaviour of 
the amplitudes is completely different. 

We should note that close to the threshold of the unitarity cut at s — 8.29, numerical errors increase and the sum 
of the LF amplitudes differs from the covariant amplitude. We have checked that the covariant amplitudes agree with 
the LF ones to at least three decimal places for s < 5.9 if we use 20 Gauss points for each dimension in the numerical 
integrals. If we increase the number of abscissas to 50, the accuracy increases to at least 5 decimal places. Thus, we 
conclude that the differences are purely due to numerical noise. 

We have also calculated the s-dependences of the contributions of the four diagrams to the form factors. In Fig. 1131 
we show the results for a fixed angle 8 — ir/2. (We do not show the results for s > 4.5, because they show little 
structure.) We have summed the contributions per diagram to each form factor and compared these sums with 
the covariant form factors. We found that the sums and the covariant results are the same, meaning that the LF 
calculation is equivalent to the covariant one. 

The contributions to the form factors are of the same order of magnitude as the sum at high s and 8. The LF 
form factors diverge for small s and 8, but their sum does not diverge since the four diagrams compensate each other, 
rendering the sum finite. The divergence at small s and 8 of the four diagrams are to some extent artefacts of the 
extraction procedure. However, they become important in a situation where e.g. the stretched box would be dropped. 
They show in a dramatic way the need to include all Fock sectors that contribute to a certain order in covariant 
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FIG. 8: The s- and ^-dependence of the form factors. Covariant calculation. 
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FIG. 9: The s-dependence of the contributions of the four diagrams to Tn in forward kinematics in 1+1 dimensions. Light-front 
calculation. 

perturbation theory. 

We have obtained our results for the form factors in LFD using the extraction procedure in Sec. IIII CI In this 
extraction procedure, we used the elements of the first row of T to obtain our form factors. Another choice for 
the four independent matrix elements of T gives different results for the contributions of the four diagrams. The 
sum, however, does not change. Furthermore, the low s- and ^-behaviour will also not change in the sum, as this 
behaviour in the individual LFD contribution is just an artefact of our extraction procedure. If one would not take 
the four-particle intermediate-state into account, in other words neglect diagram 3, the s- and ^-behaviour would be 
completely different, most dramatically at small s and 9. 

C. Light- front calculation using Pauli-Villars regularization 

In this subsection, we give our results of the LF calculation using Pauli-Villars regularization. In Figs. 1141 and 1151 
we present the results in 1+1 dimensions. All the contributions of the diagrams are convergent, and so is the sum. 
This has to be the case, since no regularization is needed in 1+1 dimensions. Furthermore, this sum is in both cases 
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— diagram 1 




diagram 2 




-- diagram 3 




— diagram 4 




— covariant 









4 5 6 7 

S 



FIG. 10: The s-dependence of the contributions of the four diagrams to Tu in backward kinematics in 1+1 dimensions. 
Light-front calculation. 




FIG. 11: The s-dependence of the contributions of the four diagrams to Tu in forward kinematics in 3+1 dimensions. Light-front 
calculation using DR2. 



equal to the covariant result. 

The results of our calculation of the Tu matrix element in forward kinematics in 3+1 dimensions at s = 5 are 
shown in Fig. 1161 The logarithmic dependence is clearly visible. The divergent parts can be written in the form 

C d lnA 2 //i 2 , (61) 

where we use the same coefficients Cd, given by Eq. H58|l . as in the case of DR 2 . Note, however, that this choice is 
not unique, as we could also subtract in addition a term of the form Cd X constant, where constant is independent 
of the diagram label d. This would not change the sum of the LF amplitudes, because the sum of the coefficients Cd 

0.07 




FIG. 12: The s-dependence of the contributions of the four diagrams to 7n in backward kinematics in 3+1 dimensions. 
Light- front calculation using DR2. 
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FIG. 13: The s-dependence of the contributions of the four diagrams to the form factors for a value 9 — n/2. Light- front 
calculation using DR2. 



vanishes. The cancellation of the logarithmic terms is clearly visible. Finally, we note that the sum is equal to the 
covariant result and that neglecting the stretched box changes the results considerably, as it is needed to obtain finite 
results. 

In Fig El we show the same results, but now we have subtracted the terms CrflnA 2 //x 2 . Again, one sees that all 
diagrams are now convergent so that this divergence is really the logarithmic divergence found by Van Iersel. 



VII. CONCLUSIONS 



We have calculated the box diagram in generalized Yukawa theory, both in the manifestly covariant formalism and 
LF quantization. In the particular case we considered, this diagram can be expressed in terms of four invariant form 
factors. 

If the calculations are carried out in 1 + 1 dimensions, both approaches lead to finite integrals and the LF amplitudes 
add up to the covariant ones. In 3 + 1 dimensions, we recovered the well-known divergences of the LF amplitudes. 
Using two completely different regularization methods, we succeeded to remove them and show that the finite parts 
of the LF amplitudes again add up to the covariant ones, while the divergent parts cancel out. There is no scheme 
dependence of the sum of the regularized amplitudes. One may change both regularization schemes in such a way that 
the individual regulated LF amplitudes in different schemes are different, but the allowed changes must not change 
their sum. There is no anomaly in the box diagram, opposed to the situation in the triangle diagram of the gauge 
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FIG. 14: The A-dependence of the four contributions to T\\ in forward kinematics in 1+1 dimensions at s = 5. Light-front 
calculation with PV subtraction. 
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FIG. 15: The A-dependence of the four contributions to 7ii in backward kinematics in 1+1 dimensions at s = 5. Light-front 
calculation with PV subtraction. 



theory. 

A zero mode occurs, that is revealed only in special kinematical situations, because it is seen only if there is no 
plus-momentum transfer. In the kinematics we adopted, it shows up for special choices of the scattering angle or the 
energy, but one might choose the scattering plane to be perpendicular to the z-axis, in which case the zero mode will 
occur for all values of the energy and the scattering angle. 

The mechanism for the occurrence of the LF singularities is easily understood. They are connected to the form 
of the dispersion relation for the LF energy. The appearance of the square of the transverse momentum in the 
numerator of the LF energy upsets the usual power counting and opens the possibility for ulraviolet divergences in 
LFD that are not present in the corresponding covariant calculation. From our results, one understands immediately 
that logarithmic divergences like the ones we have encountered in the box diagram will be present in all orders of 
perturbation theory. 

The taming of the LF singularities can be done only by including all relevant Fock sectors. In the box diagram, it 
means that one may not discard the stretched box. This does not bode well for nonperturbative methods that rely 
on rcgularization, for instance of the kernel of the bound-state equation, order- by-order. It might turn out that an 
infinite number of counterterms is needed to regulate LF Yukawa theory. 

For further study, we think that it may be useful to investigate the double box in order to learn what new singularities 
may occur. The double box would correspond to one more iteration of the covariant ladder kernel. From this study, we 
would find out whether the removal of the singularities of the box diagram is sufficient to regulate the nonperturbative 
calculations. 

We may also try to attack this problem from a completely different angle. Namely, we may just do the nonpertur- 
bative calculation for a finite value of the Pauli-Villars mass A and make aA^oo extrapolation of the observables, 
e.g. the masses of the bound states. We are convinced that this program would work if the limit of A — > oo and the 
limit of taking an infinite number of diagrams into account commute. 
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FIG. 16: The A-dependence of the four contributions to T\\ in forward kinematics in 3+1 dimensions at s = 5. Light-front 
calculation with PV subtraction. 



Since the positivity of the LF energy and longitudinal momentum is a great advantage of the LF formulation, it 
would be a disappointment if none of the cures that we would like to try would work. 



APPENDIX A: SPINORS 



We shall use in this paper helicity spinors. The reason is that in LFD the LF boosts, which are combinations 
of pure boosts and rotations, form a subgroup of the Lorent z gr oup. We get the Kogut and Soper spinors fl]| by 
boosting a rest-frame spinor using a LF boost (see, e.g., Ref. We write them down explicitly for completeness. 

(The symbols ± refer to s z = ±1/2, respectively.) 



mlf(p;+) 



p 

n + 



>/4V2p+M I V 2 P ~ m 
P r 
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«lf(p;-) 
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P l 



J 



(Al) 



The matrix elements now have the structure 
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(A2) 



The explicit forms of the matrix elements are 
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In the kinematics we adopt, we find explicitly for particle 1 

j (1) = 2£+p(l + cos6 l ) I b (l)- (E + p)psin( 



2^(E +p)(E+p cos 2My/(E +p)(E +p cos 6) 

r „ = (£+p)(£+pcos6>) +m 2 ^ _ psin<9 



2M v /( J B+p)( J B+pcos6»)' 2^(£+p)(£+pcos60' 
T x^ = (E + p)p sin g r£(l)- P(l-cosff) 



2M v / ( J E+p)(£;+pcos(9)' 2 v /(£+p)(£+pcos(9) : 
(£+p)psin6> T>V(t\ — a p(l-cos< 



r 2(!) = * , 7= == , r "(!) 



2M^(£; + p)(£; + pcos6l)' 6 2 v /(S+p)(S+ 3 9Cos6')' 
rz(1) = (E+p)(E+pcos6)~M 2 r psing 

2M v /( J B+p)(£;+pcos6l) ' 2 v /(S+39)(S+pcos6')' 

For particle 2 the matrix elements are found by substituting — p for p in the expressions for particle 1. 

A simple check on these expressions is to go to the forward limit, 6 = 0, where the matrix elements of the identity 
must reduce to 5\>\ and of the 7-matrices to p^ 6\'\/m. 



APPENDIX B: INTEGRALS NEEDED IN DIMENSIONAL REGULARIZATION 

We need explicit expressions for the integrals over k 1 - since some of the integrals are divergent. In this appendix, 
we show how to regularize these divergences using DR 2 . 

In LFD, regularization is needed in the transverse directions only. In dimensional regularization, we calculate the 
integrals as an analytic function of the dimensionality of these directions. We change the number of dimensions from 
2 to D. The integrals are changed according to 



(jfe- 1 ) +c 2 



Now, we take D = 2 — 2e to obtain 



(4tt) b / 2 T(D/2)T(s) \C 2 

For the divergent integral if, we find 
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3 4tt 
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The other integrals we need are completely regular; 
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J » - hh (B4 > 



[1] P.A.M. Dirac, Rev. Mod. Phys. 21, 392 (1949). 

[2] S.J. Brodsky, H.-C. Pauli, and S. Pinsky, Phys. Rep. 301, 299 (1998). 

[3] C.-R. Ji and C. Mitchell, Phys. Rev. D 64, 085013 (2001). 

[4] C.-R. Ji, G.-H. Kim, and D.-P. Min, Phys. Rev. D 64, 025009 (2001); Phys. Rev. D 58, 105020 (1998). 

[5] S.J. Brodsky, C.-R. Ji, and M. Sawicki, Phys. Rev. D 32, 1530 (1985). 

[6] J.H.O. Sales, T. Frederico, B.V. Carlson, and P.U. Sauer, Phys. Rev. C 61, 044003 (2000). 



21 



[7] J.H.O. Sales, T. Frederico, B.V. Carlson, and P.U. Sauer, Phys. Rev. C 63, 064003 (2001). 
[8] R.J. Perry, A. Harindranath, and K.G. Wilson, Phys. Rev. Lett. 65, 2959 (1990). 

[9] S. Glazek, A. Harindranath, S. Pinsky, J. Shigemitsu, and K.G. Wilson, Phys. Rev. D 47, 1599 (1993). 
[10] B.L.G. Bakker, M.A. DeWitt, C.-R. Ji, and Y. Mishchenko, Phys. Rev. D 72, 076005 (2005). 
[11] B.L.G. Bakker and C.-R. Ji, Phys. Rev. D 71, 053005 (2005). 
[12] M. van Iersel, Few-Body Systems 36, 133 (2005). 

[13] R. Karplus, CM. Sommerfield, and E.H. Wichmann, Phys. Rev. 114, 376 (1959) 

[14] R.J. Eden, P.V. Landshoff, D.I. Olive, and J.C. Polkinghorne, The analytic S-matrix, (Cambridge UP, Cambridge, 1966) 
[15] J.B. Kogut and D.E. Soper, Phys. Rev. D 1, 2901 (1970). 
[16] B.L.G. Bakker and C.-R. Ji, Phys. Rev. D 65, 073002 (2002). 

[17] E. Leader, Spin in Particle Physics, (Cambridge University Press, Cambridge, 2001). 
[18] K. Erkelenz, Phys. Rept. 13, 191 (1974). 

[19] N.E. Ligterink and B.L.G. Bakker, Phys. Rev. D 52, 5954 (1995). 

[20] S.J. Brodsky, J.R. Hiller, and G. McCartor, Ann. Phys. 321, 1240 (2006) 

[21] M. van Iersel, Aspects of Bound-State Calculations in Light-Front Dynamics, PhD-thesis, Vrije Universiteit, Amsterdam 
(2004) 



22 



0.06 
0.05 


— diagram 1 




: diagram 2 


0.04 


: — diagram 3 




: — diagram 4 


n m 

u.uo 


LUvdl Idl I L 


0.02 




0.01 




0.00 





2 4 6 8 10 

A 



FIG. 17: The A-dependence of the four contributions to Tu in forward kinematics in 3+1 dimensions at s = 5. In this figure 
the In A 2 //n 2 -term is subtracted. 



